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Abstract Fourier phases contain a vast amount of infor- 
mation about structure in direct space, that most statisti- 
cal tools never tap into. We address ALMA's ability to de- 
tect and recover this information, using the probability dis- 
tribution function (PDF) of phase increments, and the re- 
lated concepts of phase entropy and phase structure quan- 
tity. We show that ALMA, with its high dynamical range, 
is definitely needed to achieve significant detection of phase 
structure, and that it will do so even in the presence of a 
fair amount of atmospheric phase noise. We also show that 
ALMA should be able to recover the actual "amount" of 
phase structure in the noise-free case, if multiple configura- 
tions are used. 

Keywords Instrumentation: interferometers • Methods: 
statistical ■ Methods: numerical • ISM: structure 



1 Introduction 

Observations of the interstellar medium (ISM) reveal highly 
complex, fractal-like structures [8]|7]. The self-similar hier- 
archy of these structures, over four decades, is thought to 
spring from the interplay of turbulent motions Ifl2ll and self- 
gravitation 0. To understand this interplay, one therefore 
needs a quantitative description of the observed structures. 
Most of the statistical tools used to this end are more or less 
derived from the power spectrum [6], which is given by the 
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squared amplitudes of Fourier components. Yet, a simple 
numerical experiment performed by [5] shows that essential 
structural information lies in the Fourier-spatial distribution 
of the phases. 

In the following, we present some of the notions used 
to exploit this information (section 0, and their practical 
implementation (section [3}. We then consider the ability of 
ALMA and other arrays to detect and measure phase struc- 
ture information in real time (section @}. We conclude by 
giving some future perspectives (section[5]l. 



2 Fourier phase analysis 

The importance of Fourier phases in terms of structure has 
been recognized by various studies [ 15, 13,5 1. Since the in- 
formation sought lies in the Fourier spatial distribution of 
phases, Scherrer et al. 03 1 suggested considering the statis- 
tics of phase increments As<p{k) - <f>(k + 6) - 4>(k) between 
points separated by a given lag vector 6 in Fourier space. 

In a field for which Fourier phases are uncorrelated, such 
as fractional Brownian motions (fBm£] fl6l . phase incre- 
ments are uniformly distributed over [—re, n], for any lag vec- 
tor 6. At the other end of the spectrum is the case of a single 
point source, for which the PDF of phase increments is a 
delta function. In between those extremes, the PDF of phase 
increments presents a single wavelike oscillation (See Fig.Q] 
for an example), which may be seen as a signature of phase 
structure. 

A quantitative measure of the distribution's departure 
from uniformity is phase entropy [13], 
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1 These are random fields characterized by a power-law power spec- 
trum and random phases. 



2 



which reaches its maximum value Sq = ln(2w) for ffims. 
It is therefore convenient to consider the positive quantity 
Q(6) = So - S(6), which we dub phase structure quantity, 
and which may be directly computed on the histograms of 
phase increments. 
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Fig. 1 Top : Column density of a 512 3 weakly compressible hydrody- 
namical turbulence simulation obtained by Porter et al. 1141 . used here 
as a model brightness distribution for phase structure analysis. Bottom 
: Histogram of phase increments for this field, with S = e x (unit vec- 
tor along the k x axis in Fourier space) and n = 50. The dotted line 
represents the uniform distribution. 



3 Phase structure quantity in practice 

For a finite-sized image, histograms of phase increments do 
not perfectly sample the underlying PDFs. Phase structure 
quantities Q associated with these distributions should be 
distinguished from those Q found by numerical integration 
of the histogram^, which depend on the number p of avail- 

2 To give an idea, for the histogram shown on Fig. [T] we have 
(2(e. r ) = 9.8xlCT 3 . 



able increments and the number n of bins. The difference 
between the two can become significant for Q <sc 1. 

In particular, finding Q ^ does not guarantee that phase 
structure is present in a given field (Q + 0), so that the de- 
tectability of phase structure depends on the threshold of Q 
above which there is a given probability (say 0.99) that an 
image deviates significantly from a "structureless" field. 

The procedure is described in [11] and is largely based 
on results from [4 |. In short, the end result is that the thresh- 
old of Q depends on n and p, and may be found using well- 
known^- 2 statistics. 

The influence of n and p on the reliability of Q may also 
be studied numerically, using fractional Brownian motions. 
Unsurprisingly, while Q — for these, Q increases as the 
size of the image decreases, and as the number of bins in- 
creases. 



4 Application to interferometric observations 

In the ideal case, interferometers sample the Fourier trans- 
form of observed brightness distributions, and allow direct 
measurement of phase increments. Since this can be done 
as the Earth rotates, we may look for the minimum observ- 
ing time required to detect a significant phase structure in 
the data. To focus on the problem of statistical estimation 
described in the previous section, we shall not consider pri- 
mary beam attenuation nor regridding issues. These simpli- 
fications are discussed in more detail in IfTTII . 

To estimate the ability of ALMA to detect and measure 
phase structure, we proceed as follows: A model brightness 
distribution is taken as input to a simple interferometer sim- 
ulator, which is based on the characteristics of ALMA and 
uses the array configurations optimized by Boone [ 1 1. The 
instrument tracks the source as long as it remains above a 
minimum elevation of 10°. The output maps, for which no 
deconvolution is performed, yield values of Q as a function 
of integration time, with 6 and n fixed. 

The model brightness distributions used are the one of 
Fig. Q] and a field with the same power spectrum, but with 
random phases. For comparison, we have also considered 
configurations taken from current arrays, such as the Plateau 
de Bure (PdB) and the VLA, fictitiously located at the same 
geographical coordinates as ALMA, and observing the same 
source. 

As the observation is carried out, more and more Fourier 
phases are measured and p increases. The question is whether 
this allows to bring down the upper limit discussed in sec- 
tion|3] below the measured phase structure quantities, to en- 
sure positive detection. The results are summarized on Fig- 
ures [2] to H] which show the evolution of Q(e x ) as a func- 
tion of integration time. Fig. [2] shows that the number of 
phase increments measured by the Plateau de Bure in its 



3 



1.0 



% 0.4 

■1; 

en 
o 

* 0.2 



0.0 

2 4 6 8 10 

htegrotion time (h) 

Fig. 2 Evolution of measured Q(e x ) with integration time, for the B 
configuration of the Plateau de Bure. The black solid line corresponds 
to the turbulent brightness distribution, and the grey solid line to the 
random-phase brightness distribution. The dotted line represents Q(e x ) 
for the complete turbulent brightness distribution, and the dashed line 
represents the evolution of the theoretical upper limit (lying above the 
plotted range here). 

B configuration is insufficient to detect phase structure, as 
the curves for turbulent and random-phase brightness dis- 
tributions are indistinguishable from one another. The same 
conclusion prevails for other configurations of this instru- 
ment and other lag vectors. On the contrary, Fig. [3] shows 
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Fig. 3 Same as Fig.|2] but for the D configuration of the VLA. 

that the VLA allows such a detection, since the measured Q 
becomes larger than the theoretical upper limit, after about 
6 hours of integration. Long before that, however, we get 
a hint that phase structure is present in the field, since the 
curves for both model brightness distributions go apart after 
less than twenty minutes. This diagnosis can be performed 
in real time by drawing random phases for the visibilities as 
they are measured. ALMA gives even better results (Fig.@}. 
In its E configuration and in our case, a short integration 
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Fig. 4 Same as Fig. [2] but for the E configuration of ALMA. 

time of about twenty minutes is enough to conclude on the 
presence of phase structure. However, the final value of Q 
obtained is not equal to the phase structure quantity mea- 
sured on the model brightness distribution. This is due to 
the fact that only 24% of the 512 x 512 Fourier phases are 
measured by this configuration. 

Using more extended configurations, one should be able 
to measure the Fourier components lying outside the radius 
covered by the E configuration, and therefore hope to re- 
cover the correct value of the phase structure quantity by 
combining visibilities from multiple configurations. Fig. [5] 
shows the evolution of the measured Q(e x ) with integration 
time, using this approach)! It appears that the Fourier plane 
coverage achieved by ALMA will allow measurement of the 
actual value of the phase structure quantity for the observed 
field, while the VLA fails. 

Finally, to assess whether atmospheric phase noise would 
prevent detection of phase structure, we introduced a mask 
giving the refractivity field above the instrument. We as- 
sumed this mask to be a 200-m thick layer of frozen Kol- 
mogorov turbulence being transported along the east-west 
direction at 2 m.s -1 , and normalized it so that the rms phase 
noise <x for a pair of antennae observing the zenith and sep- 
arated by a baseline d = 100 m should be one of a few 
specific values, namely 15°, 45° and 90°. According to 
and using the scaling relation given by [9|, noise levels at 
Chajnantor vary typically from <x ~ 14° to ctq ~ 57°. 

Integration of the refractivity field along the different 
lines of sight for each antenna as the observation is per- 
formed yields phase delays, which are then correlated to ob- 
tain the atmospheric phase noise for each pair of antennae, at 
all times. Fig. [6] shows the evolution of the measured Q(e x ) 
for ALMA in its E configuration. 



3 The integration time r is to be understood per configuration, and 
the total time of integration is ^configurations x r. 
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Fig. 5 Evolution of measured Q(e x ) with integration time for an ob- 
servation using all configurations of the instrument in turn. The black 
solid line corresponds to the six configurations of ALMA, and the grey 
solid line to the four configurations of the VLA. The dotted line repre- 
sents the value of Q(e x ) for the whole field. 
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Fig. 6 Evolution of measured Q(e x ) with integration time in the pres- 
ence of atmospheric phase noise (solid lines, with o"o specified next to 
each curve). The array used is the E configuration of ALMA. The dot- 
ted line represents Q(e x ) for the whole field, and the dashed line shows 
the theoretical upper limit. 



It appears that in this case, the presence of phase struc- 
ture can be easily detected in the presence of a fair amount 
of atmospheric phase noise. Indeed, even a rms phase fluctu- 
ation of (To = 90° is insufficient to bring the measured phase 
structure quantity below the upper limit. Consequently, phase 
structure will undoubtedly be detected by ALMA without 
any phase correction, although the use of dedicated water 
vapor radiometers, as is planned, should allow for an effec- 
tive decrease of the atmospheric phase noise by a substantial 
factor [ 10 1, making it possible to actually measure the phase 
structure quantity for the observed field. 



5 Perspectives 

In the context of interferometry, a more elaborate use of 
phase information would be to keep track of the phase mea- 
sured by each baseline as a function of time, and to compute 
phase increments along the baseline's track. This should re- 
duce contamination by atmospheric phase noise, but would 
require a shift in the phase structure information formalism, 
since, in this approach, the lag vector 6 is no longer a control 
parameter, but a function of time and of the baseline. 

Another possible extension of this work is the inclu- 
sion of the kinematic dimension, which is accessible through 
ALMA's high spectral resolution receivers. It may well be 
that phase analysis applied to individual channel maps should 
prove a valuable tool for assessing the three-dimensional 
structure of velocity fields. 
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